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, Late-type giant stars in the evolutionary stage of the asymptotic giant branch increasingly lose mass via comparatively slow but 

1/^ ■ dense stellar winds. Not only do these evolved red giants contribute in this way to the enrichment of the surrounding interstellar 

medium, but the outflows also have a substantial influence on the spectro-photometric appearance of such objects. In the case of 
carbon-rich atmospheric chemistries, the developing cool circumstellar envelopes contain dust grains mainly composed of amorphous 
carbon. With increasing mass-loss rates, this leads to more and more pronounced circumstellar reddening. With the help of model 
. calculations we aim at reproducing the observational photometric findings for a large sample of well-characterised galactic C-type 

' Mira variables losing mass at different rates. We used dynamic model atmospheres, describing the outer layers of C-rich Miras, 

which are severly affected by dynamic effects. Based on the resulting structures and under the assumptions of chemical equilibrium 
, as well as LTE, we computed synthetic spectra and synthetic broad-band photometry (Johnson-Cousins-Glass BVRIJHKL'M). A set 

of five representative models with different stellar parameters describes a sequence from less to more evolved objects with steadily 
Q ' increasing mass-loss rates. This allowed us to study the significant influence of circumstellar dust on the spectral energy distributions 

< I ' and the (amplitudes of) lightcurves in different filters. We tested the photometric properties (mean NIR magnitudes, colours, and 

, amplitudes) and other characteristics of the models (mass-loss rates, periods, and bolometric corrections) by comparing these with 

^ ■ the coiTesponding observational data adopted from the literature. Using different kinds of diagrams we illustrate where the models are 

located in a supposed evolutionary sequence defined by observed C-type Mira samples. Based on comparisons of galactic targets with 
empirical relations derived for C stars in the Large Magellanic Cloud we discuss the relevance of metallicity and excess carbon (C~0) 
for the development of dust-driven winds. Having investigated the dynamic model atmospheres from different (mainly photometric) 
^ . perspectives, we conclude that our modelling approach (meaning the combination of numerical method and a suitable choice of model 

in ' parameters) is able to describe C-rich long-period variables over a wide range of mass-loss rates, i.e., from moderately pulsating 

objects without any dusty wind to highly dust-enshrouded Carbon Miras. Thus, we can trace the observed sequence of C-type Miras, 
, which is mainly determined by the mass loss. 

. 
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1. Introduction objects. In Paper II we presented one selected dynamic model 

atmosphere that is representative of a typical C-type Mira vari- 
This work IS the third m a series of publications dealing with ^^j^ j^^j^g ^^^^ ^„ intermediate rate and that accounts for 
synthetic spectra and photometiyof carbon-nch late-type giants, ^j^^^^ dynamic effects. Simulating the pulsation with a vari- 
H In the first part, Annger et al. (|2009] below Paper I) presented a ^^le inner boundary condition results in large-amplitude photo- 
^y. grid of hydrostatic, dust-free model atmospheres and the cor- ^^^^^^ variations similar to observed lightcurves. The circum- 
responding photometnc results. It was found that the synthetic ^^^jj^j. ^^^^ composed of amorphous carbon (amC) grains leads 
photometry fits observational data reasonably well for objects even higher amplitudes, as well as to redder and more realis- 
with eff^ective temperatures above ^2800K. colours compared with the hydrostatic case. Both findings 
In the second part, Nowotny et al. (12011a) from now on are characteristic of long-period variables (LPVs), which are 
Paper II) extended this study to more evolved C-rich stars for evolved red giant stars in the evolutionary phase of the asymp- 
which dynamic effects become important. With outer atmo- totic giant branch (AGB). We investigated in detail the resulting 
spheric layers strongly influenced by (i) shock waves caused synthetic lightcurves and variations in diff-erent colour indices as 
by pulsations of the stellar interiors and (ii) the development of ^ function of phase 0boi • The comparison proved that these mod- 
dusty stellar winds, a diff-erent modelling approach was needed glling results resemble the observed fight variations of the C-rich 
when trying to reproduce the observational properties of such Mira RU Vir reasonably well 



* Appendices are only available in electronic form at The choice of model in Paper II, i.e. one with an intermedi- 

http : //www . aanda . org ate mass-loss rate, was motivated by two reasons: (i) to study the 
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effects of pulsation and dust in a typical Mira, and (ii) to show 
that our modelling approach gives realistic results for stars where 
the emergent flux is a mixture of photospheric and circumstellar 
contributions. Extreme cases, i.e. dust-free, only mildly pulsat- 
ing stars on the one hand, or objects with optically thick dusty 
envelopes on the other, can be reproduced rather well with hy- 
drostatic atmosphere models (Paper I) or earlier generations of 
dust-driven wind models (e.g. Hofner & Dorfi 1997, Winters et 
al.|2000). However, to achieve a realistic representation of both 
the stellar atmosphere and the dusty outflow, it is necessary to 
apply models that combine time-dependent dynamics, detailed 
non-equilibrium dust formation, and frequency-dependent radia- 
tive transfer (Hofner et al. .2003. Nowotny et al. 2010. Paper II). 

In the present paper we study the influence of stellar parame- 
ters on dust-driven outflows and the resulting photometric prop- 
erties, illustrated by a sequence of models representing different 
phases of AGB evolution. We demonstrate that our modelling 
approach allows us to cover a wide range of observed properties 
by adjusting a few well-defined input parameters. In this way 
we can compute realistic synthetic photometry for models rang- 
ing from dust-free pulsating atmospheres through objects with 
intermediate mass-loss rates to stars with optically thick dusty 
outflows. 

For clarity, it should be mentioned here that the emergent 
wind properties (dust characteristics, mass-loss rates, velocities, 
optical depths, etc.) are results - not additional input parameters 
- of our models that link stellar atmospheres and circumstellar 
envelopes in a self-consistent way. In contrast to other often-used 
approaches (e.g. DUSTY code) where the outcome of a dust ra- 
diative transfer calculation is controlled by the choice of input 
parameters (i.e. properties of the circumstellar dust shells) quite 
directly, the relation between basic stellar properties and syn- 
thetic observables is non-trivial in our ab initio modelling. The 
non-linearities inherent to the self-consistent dynamical mod- 
elling make a direct fitting of specific objects by fine-tuning of 
input parameters complicated. Notably, it is not possible to cre- 
ate arbitrary combinations of atmospheres and circumstellar en- 
velopes, and certain sets of stellar parameters will not lead to 
outflows at all. On the positive side, however, our models have a 
significantly higher predictive power since observable properties 
that are linked in real objects cannot be adjusted independently 
in the models. 

In this paper we restrict the discussion of trends of photomet- 
ric properties with fundamental stellar and pulsation parameters 
to a small group of models, in order to present them in some de- 
tail. A paper covering photometric results based on a significant 
fraction of the grid of dynamic model atmospheres presented by 
Mattsson et al. (2010) is currently in preparation. The aim is to 
make the models available for interpreting of individual stars, as 
well as for use in stellar evolution and stellar population mod- 
elling, by providing a consistent set of wind properties, together 
with synthetic spectra and photometry for a wide range of pa- 
rameters. 

The model characteristics and the computational details are 
outlined in Sect.|2] while in Sect.[3](as well as App.lAl) the com- 
parative data compiled from the literature is presented. In Sect.|4] 
we describe how the spectral energy distributions (SEDs) and 
the photometric variations are influenced by the model param- 
eters, especially by the differing properties of the circumstellar 
dust envelopes. In Sect.|5]a comparison of various quantities with 
the corresponding measurements for a sample of galactic Carbon 
Miras as found in the literature is carried out. Thereafter, we 
study in Sect.|6]how the models fit into an evolved stellar popu- 
lation (also containing among various red giant stars the targets 



of our atmospheric modelling, namely C-rich LPVs with mass 
loss), and discuss the importance of free carbon (C-O) for the 
dust formation. 



2. Model atmospheres and radiative transfer 

2.1. Combined atmosphere and wind models 

For the study presented here we used the dynamic model atmo- 
spheres of Hofner et al. (2003), which cover the two major dy- 
namic aspects relevant to the outer layers of evolved red giants, 
namely the effects of pulsation and mass loss via a dust-driven 
wind. To this end, the equations for hydrodynamics, frequency- 
dependent radiative transfer and dust formation are solved si- 
multaneously. This modelling approach provides a realistic and 
self-consistent description of the dynamic stellar atmosphere as 
well as the regions of the outflow, which is triggered by radiation 
pressure on newly formed amC grains. Such atmospheric models 
are well suited to simulate mass-losing LPVs with carbon-rich 
atmospheric chemistry (C/O > 1). 

For more information about the numerical and physical de- 
tails we refer to a number of previous publications dealing 
with this (e.g. Hofner et al. 2003; Gautschy-Loidl et al. l2004t 
Nowotny et al. 2005a. 2005b, 2010, 2011ai Mattsson et al. 
2010). Graphical representations of typical models can be found 
in Nowotny et al. (120101 1201 lal) . These include plots of the tem- 
porally varying radial atmospheric structures (p, T, p, u, fc) at 
different phases 0boi during the pulsation cycle and illustrations 
of the movements of atmospheric layers at different depths. 

In Paper II we selected one specific dynamic model atmo- 
sphere (named model S) and studied in detail the influence of 
circumstellar dust on the photometric appearance as well as the 
photometric variations throughout a pulsation cycle. This model, 
also included here (cf. Table[T]i, resembles a typical carbon-rich 
Mira with an intermediate mass-loss rate M of approximately 
a few 10"'' Mq yr"'. In addition, we investigate a larger sample 
of models. The parameters of the respective hydrostatic initial 
models are listed together with the resultant wind properties in 
TablelU 

Although existing objects served as rough reference points 
for the chosen parameter combinations, the five models should 
not be regarded as a dedicated fit (for a discussion see Nowotny 
et al. I2005bl) for any of the targets used for the comparison 
in Sect.|5] Instead, the aim was to retrace with the help of 
characteristic models the sequence of galactic C-type Miras 
recognisable in observational studies (see below). Therefore, the 
model parameters were varied to proceed from less to more 
evolved objects. This means that while the effective temperatures 
Ti, and the surface gravities decrease for the model series 
T^R^S-^Cl — > C2, increasing values were chosen for the 
luminosities L*, the radii the C/O ratios, and the pulsation 
periods P. Since one of the significant effects along the observed 
sequence is the increasingly relevant circumstellar reddeningQ 
we also applied models characterised by increasing mass-loss 
rates M. The latter are dependent on the piston velocity ampli- 
tude AMp (see Nowotny et al. 1201 Oi l. With appropriately chosen 
values for these input parameters we can influence the arising 
outflow to some degree. However, occurrence and intensity of 



' Throughout the paper the term "reddening" is used in two different 
senses: (i) the colour of a star changing due to the dust grains in the cir- 
cumstellar envelope surrounding it, and (ii) the colour becoming redder 
because of the interstellar extinction of the light on its way from the 
object to us. 
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Table 1. Characteristics of the dynamic model atmospheres for pulsating and mass-losing C-rich AGB stars used for the photometric 
modelhng. 



Model: 


T 


R 


S 


CI 


C2 


(1) 


T 

1^* 




5200 


7000 


10 000 


10000 


10000 




A/I 
™* 


L^oJ 


1 


1 


1 


1.5 


1.5 




T^, 


[K] 


3000 


2800 


2600 


2400 


2400 




[Fe/H] 


[dex] 



















C/O 


by number 


1.1 


1.4 


1.4 


1.69 


2.38 


W 




rp 1 


268 


355 


493 


580 


580 






[AU] 


1.24 


1.65 


2.29 


2.70 


2.70 




log (g* [cms-2]) 




-0.42 


-0.66 


-0.94 


-0.91 


-0.91 


(Hi; 


p 




295 


390 


490 


640 


640 




P 


[km s-i] 


4 


5 


4 


6 


6 




/l 




1 


1 


2 


2 


2 




AnJbol 


[mag] 


0.41 


0.53 


0.86 


1.49 


1.49 


(iv) 


<M) 


[Moyr-'] 




3.5 X 10-'' 


4.3 X 10-« 


9.1 X lO-'^ 


1.4x10-5 




(u) 


[km s-'] 




15 


15 


17 


28.6 




(fc) 






0.33 


0.28 


0.43 


0.58 




(Pd/Pg) 


[10-3] 




0.75 


0.63 


1.22 


3.23 


(V) 


Tdusi (0.55 yum) 






3.26 


4.00 


11.06 


19.45 




Tdu,sl(l-6yUm) 






1.01 


1.24 


3.42 


6.02 




Tdiist (3.8yum) 






0.37 


0.45 


1.25 


2.20 



Notes. Listed are (i) parameters of the hydrostatic initial models, (ii) quantities derivable from these parameters, (iii) attributes of the inner 
boundary conditions (piston) used to simulate the pulsating stellar interiors as well as the resulting bolometric amplitudes Amtoi, (iv) properties of 
the resulting wind, and (v) dust optical depths of the formed circumstellar envelope. 

The notation follows previous papers (Hofner et al. 120031 Nowotny et al. 120101 1201 la ll: i?^. - stellar radius of the hydrostatic initial model 
calculated from the luminosity and temperature Ti, via the relation =AnR\. crTl; P, Ai(p - period and velocity amplitude of the piston at 
the inner boundary; - free parameter to adjust the luminosity amplitude at the inner boundary; (M), {ii} - mean mass-loss rate and outflow 
velocity at the outer boundary; ( f^) - mean degree of condensation of the element carbon into amC dust grains at the outer boundary; (pd/Pg) 
- mean dust-to-gas ratio at the outer boundary (cf. Eq. (2) of Hofner & Dorfi 1997); Tdust - optical depth of the dust envelope at the specified 
wavelength (V-, H-, L'-band) based on radial structures at (f)f,„i a 0.25 (see Sect. l4.2t . The radial coordinates in this work are plotted in units of the 
corresponding stellar radii of the hydrostatic initial models, calculated from their luminosities and temperatures (as given in the table) 
via the relation L^, = 4nR\ aT* . 



the wind are an outcome of the modelling. In our case, we ar- 
rived at models as different as model T (dust-free, pulsating at- 
mosphere without any outflow) and model C2 (substantial M ex- 
ceeding IQ-^ Moyr-' and leading to an optically thick circum- 
stellar dusty envelope). They presumably mark extreme values 
of the mass-loss sequence found observationally, while the pre- 
viously studied model S is located roughly in the middle. 

Model R was adopted from Hofner et al. (I2003I I. while mod- 
els T and S were taken from Gautschy-Loidl et al. (2004). The 
latter model was also used in various other studies (Nowotny et 
al. I2005al l2005bl l20T0l ISoTTal Paladini et al. ,2009) . Computed 
specifically for this study were the last two models, CI and C2. 

2.2. Synthetic spectra and photometry 

Snapshots of the atmospheric structures at various phases dur- 
ing several periods (see, for example. Fig. A. 2 in Paper II) were 
then used for the detailed a posteriori radiative transfer calcula- 
tions. To this end we closely followed the numerical approach 
described in Paper II, for the details and an extensive list of ref- 
erences we refer to this previous paper. 

On the basis of the radial structures we applied our opac- 
ity generation code COMA (Aringer 120001 Paper II) to compute 
atomic and molecular abundances for every depth point of the at- 
mospheric models with equilibrium gas chemistry routines. The 
principle element abundance pattern of solar composition was 
modified concerning carbon in two ways. First, the general C 



abundance was changed according to the C/O ratios as speci- 
fied for the different models (cf. Table[T]i. Second, the depletion 
of C in the course of amC dust grain formation was also taken 
into account. Thereafter, opacities for all atmospheric layers and 
every wavelength point were calculated assuming conditions of 
LTE. All the sources needed for an appropriate spectral synthesis 
were included: (i) continuous absorption, (ii) atomic and molec- 
ular line contributions, and (iii) opacities due to particles of amC 
dust under the assumption of the small particle limit. Again we 
followed the opacity sampling approach of Paper II to account 
for all molecular species necessary for a proper representation 
of the synthetic spectra (CO, CN, C2, CH, HCN, C2H2, C3). 

Based on the resulting arrays of combined opacities - with 
a spectral resolution of R-A/AA= 18 000 - spherical radiative 
transfer was solved. The resulting synthetic spectra in the wave- 
length range of 0.2-25 yum can directly be used to compute fil- 
ter magnitudes. For the presentation of the spectra in Fig. [Tithe 
spectral resolution was decreased to R = 360. 

Following Paper II, we chose for the synthetic photome- 
try the Johnson-Cousins UBVRI system described by Bessell 
(II99OI 1 as well as JHKLL'M from the homogenised Johnson- 
Glass system described in Bessell & Brett (I1988I I. The filter 
transmission curves given by these authors were then convolved 
with the synthetic spectra. Applying adequate photometric zero- 
points for every filter (adopted from Bessell et al. 11998) led to 
calibrated synthetic absolute magnitudes. 
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Again the Johnson-Cousins-Glass or Bessell photometric 
system was defined as the standard system for the whole work 
presented here, and the observational results (Sect.|3] App.lAli 
were transformed where required, allowing for a direct comp ar- 
ison with the modelling results. 

3. Observational comparative data 

As the main observational reference we used the compilation 
of time-series, near-infrared (NIR) photometry published by 
Whitelock et al. ('2006; below W06). They presented compre- 
hensive data sets for a large sample of galactic carbon-rich AGB 
stars. Out of the 239 LPVs listed we only adopted the subsam- 
ple with pronounced variablity, i.e. the well-characterised Miras 
in their Table 6. According to Whitelock et al. these objects are 
characterised by clearly periodic lightcurves and peak-to-peak 
K amplitude^ larger than 0.4™^, in their nomenclature this cor- 
responds to the variability class In. This group of Mira variables 
is more appropriate for an interpretation by our modelling ap- 
proach, while the remaining objects (non-Miras, variability class 
2n) were not included here. In Paper II we made a detailed com- 
parison of the lightcurves of the C-rich Mira RUVir with the 
corresponding photometric variations of model S (Table[T]- This 
is also followed up here to some extent. However, beyond that 
we made use of the complete sample of observed field C-type 
Miras. This allows for an extended and more general comparison 
with models of different parameters. To this end, we adopted the 
mean values for colour indices from Whitelock et al. rather than 
the individual measurements of the time-series, NIR photometry 
also given. The authors provide photometric results akeady cor- 
rected for interstellar reddening as well as useful (estimates for) 
other quanitities which can be compared to our models. The lat- 
ter include periods, mass-loss rates M, photometric amplitudes 
AA", bolometric magnitudes mboi, or distances (which are essen- 
tial to shift the apparent magnitudes to an absolute scale; e.g. 
Fig.fTTTl. The photometric results of Whitelock et al. (2006i were 
not transformed to our standard system (Sect. 12. 21 as the differ- 
ences between the SAAO photometric system and the Bessell 
filters are small according to Bessell & Brett (119881 see their 
Table I). 

We also adopted various other data from the literature for 
the comparisons below, a detailed description of those including 
references and the necessary post-processing can be found in 
App.|Al 

4. IVIodelling results - the influence of increasing 
mass loss 

4.1. Spectral energy distributions 

In Paper II we demonstrated with the help of model S the difi'er- 
ence between hydrostatic objects and dynamic objects (charac- 
terised by pulsating stellar interiors and winds) concerning the 
radial atmospheric structures, the resulting spectra and the pho- 
tometric appearance. We also outlined the role of different opac- 
ity sources for the spectral energy distributions, in particular we 
showed the effect of the carbon dust particles in the wind. 

Here we continue the theoretical investigations with the ex- 
tended set of dynamic model atmospheres of Table[T] These 



^ In Fig.[8]the amplitudes of Fourier-fits to the lightcurves (first-order 
sine curves) are plotted, which can be smaller than the corresponding 
peak-to-peak; amplitudes. 



models have different properties, with one of the most signifi- 
cant being the mass-loss rate. Its influence on the synthetic low- 
resolution spectra is sketched in Fig.lT] for which we used three 
models to illustrate the progression from no mass loss to highly 
dust-enshrouded. For each model we plot the spectra for a max- 
imum and a minimum phase to mark the variations during the 
pulsation cycle. 

The top panel of Fig.[T]shows the results for model T, which 
experiences no dust formation and, thus, no stellar wind oc- 
curs. The dynamical calculation leads to spectra that are rel- 
atively close to the spectrum based on the hydrostatic initial 
model (which is the starting point for the dynamic calculation; 
cf. Paper II). While all spectra are dominated by the prominent 
molecular features that are characteristic of late-type giants, the 
pulsation - simulated by the varying inner boundary - results 
in small variations of the spectrum around the hydrostatic case 
during a light cycle. In addition to the overall changes in the flux 
level because of the varying luminosity input, we find changing 
feature intensities due to the periodically modulated atmospheric 
structure (cf. Fig. 3 in Nowotny et al. l20lOi ). 

The coiTesponding results for model S are shown in the mid- 
dle panel of Fig.[T] The stellar wind containing amC grains, 
which emerges when the conditions for dust condensation are 
fulfilled, has considerable influence on the spectral appearance 
of the object. The first recognisable effect is the attenuation of 
the molecular features in the spectra based on the dynamical 
phases when compared to the hydrostatic spectrum. The second 
effect caused by the circumstellar dust is the redistribution of 
the emitted flux from the visual to the IR. This results in spec- 
tra which are clearly different from the hydrostatic one for every 
phase (/>boi during the whole pulsation cycle. 

The bottom panel of Fig.[T] displays the dramatic spectral 
change for even more evolved objects. Model C2 represents a 
star with pronounced mass loss leading to an optically thick 
dust envelope. The molecular contributions in the spectra ap- 
pear substantially weakened, except for some features (e.g. the 
characteristic one at « 3 /^m due to C2H2 and HCN), which are 
still recognisable, although they also have lower intensities. The 
whole SED is completely shifted to infrared wavelengths and 
differs significantly from the one for the corresponding initial 
model (pure photosphere, no wind). In addition, severe changes 
in flux during the light cycle can be seen. 

Figure[T] illustrates that our modelling approach allows de- 
scribing carbon-rich LPVs covering a range of quite different 
mass-loss rates. Starting with objects showing spectra charac- 
terised only by photospheric features (atomic, molecular), we 
proceed to objects for which the spectral appearance is almost 
completely determined by the circumstellar envelope resulting 
from the outflow that contains amC dust particles. This effect is 
also demonstrated in Fig.|2]where we show the SEDs of the three 
selected models of Fig.[T]using calibrated filter magnitudes com- 
puted on the basis of the corresponding spectra. Only minor de- 
viations (< l^^s) from the hydrostatic case can be found for the 
two dynamic phases of model T. However, the reddening due to 
the dust becomes severe for the other models. On the one hand, 
the objects appear fainter by several magnitudes in the visual due 
to the absorption caused by the circumstellar dust. The increase 
in this dimming towards shorter wavelengths reflects the proper- 
ties of amC dust (cf. insert of Fig.lTJ and can be larger than the 
periodic photometric variations (model C2). On the other hand, 
one can also recognise circumstellar dust emission for the filters 
L'M. The change from dimming to additional emission occurs 
at around 2 yum, which is one of the reasons why the K filter is 
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model T 



hydrostotic initiol model 

maximum phose 
minimum phose 




2 3 5 7 10 
A [/tinn] 

Fig. 1. Synthetic low-resolution spectra (R - 360) for three dy- 
namic model atmospheres as denoted in the legend. Plotted are 
spectra based on the corresponding hydrostatic initial model as 
well as two representative phases of the dynamic calculation 
(<^boi ~ 0.0 and 0.5). The insert in the middle panel shows the ab- 
sorption data for amorphous carbon dust from Rouleau & Martin 
( I1991I I as applied in the modelling. The wavelength ranges of 
the used broad-band filters FWHM of the responses given in 
Bessell 1990 and Bessell & Brett , 1988 J are marked for orienta- 
tion purposes in the lower panel. 



often used in variability studies (My^:-/" -relation) of LPVs (e.g. 
Whitelock l20T2] i. 



4.2. Optical depths of the dust envelopes 

To characterise the circumstellar envelopes we computed 
monochromatic dust optical depths Tdust at different representa- 



-15 
-10 
-5 

5 
10 
15 
20 




hydrostotic initiol model 

maximum phase 
minimum phose 



model T 




Fig. 2. Synthetic SEDs described by calibrated photometric 
magnitudes in the used broad-band filters resulting from the 
spectra of Fig. [T] 

tive wavelengths by radially integrating inwards from the outer 
boundary of the model rmax^ 



Ti.dustCO 



KA,dust(r') p(r') dr' 



(1) 



This corresponds to the central beam with impact parameter 
p-Oin the framework of spherical radiative transfer (e.g. Fig. 4 
ofYorke 1980). 

The middle panel of Fig.|3]shows how the optical depth of the 
dusty envelope varies with time for model S. Superposed on the 
more or less constant degree of dust absorption due to the steady 
outflow we also find a temporary increase in Tjust around phases 
of light minimum where enhanced dust formation takes place 
in a relatively narrow region at ~ 2 (cf. upper panel of the 
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-0.5 0.0 0.5 1.0 1.5 



-4 -, 



-2 - 



C7> 

E 



1 1 1 1 1 1 1 1 

''"•ms^ll'"' 


1 1 1 1 1 1 1 1 1 
■ • •. „oa ' ' * 

•'Si....-' 


1 1 1 1 


~_ V (no dust) 






V 












" , , , 1 , , , , 


, , , , 1 ", , , , 


, , , , " 



-0.5 



0.0 0.5 
phase 



1.0 



1.5 



Fig. 3. Illustration of the effect of dust formation on the pho- 
tometric appearance of the star, demonstrated with the help of 
model S. Upper panel: Movement of atmospheric layers of the 
model over about two light cycles, colour-coded is the degree of 
condensation of carbon into dust grains. Middle panel: Optical 
depths of the dust envelope as a function of time, calculated for 
the central wavelengths of the filters V and K following Eq. ([T]i. 
Lower panel: Resulting synthetic lightcurves for the same filters 
(filled circles), together with the corresponding photometry for 
which the dust opacities were not included (open circles). Note 
that the K-lightcurves were shifted by 5™^ to fit into the panel 
(cf. Fig. 7 in Paper II). 



is also reflected in the lightcurves shown in the lower panel of 
Fig. [3] There we plot synthetic photometry resulting from the 
full spectral synthesis (filled circles), as well as the correspond- 
ing results of a calculatioij^ where the dust opacities were not 
taken into account in the a posteriori radiative transfer (open cir- 
cles). While the photometric variations of the latter mainly re- 
flect the changing luminosity input introduced by the variable 
inner boundary of the model (with minor contributions due to 
varying molecular feature intensities), the former are in addition 
influenced by the dust effects. For the K-filter only small dif- 
ferences are recognisable at phases where a new dust shell is 
formed, whereas the V-lightcurve is severly altered by the vary- 
ing amount of dust absorption. On top of the general dimming 
by «i 4™° because of the dusty outflow (visible e.g. at ^boi = 0.0) 
an additional drop in V is clearly noticeable around 0boi ~ 0.6. In 
this way, the (time-dependent) absorption by dust grains in the 
circumstellar envelope is responsible for the generally reddened 
colours of such objects and also for the variations in colour in- 
dices, such as iV -K). 

Some cycle-to-cycle variations concerning the dust forma- 
tion and, thus, the V-lightcurve can be recognised in Fig.|3] This 
effect is relatively moderate for model S but can be much more 
pronounced for other parameter combinations (see App.lCi. 

Although we have shown that the dust optical depth is a 
time-dependent quantity we only give one characteristic num- 
ber in Table[r| The values listed there correspond to phases of 
0boi ~ 0.25 to avoid the intervals of extreme dust formation dis- 
cussed above. In this way, the Tdust should be regarded as lower 
limits and serve as representative quantites to compare the in- 
dividual models. We note (i) the monotonic increase of optical 
depths for the model sequence R— »S— >C1^C2 (i.e. increas- 
ing mass-loss rates), and (ii) the general decrease of Tdust when 
going from the visual to the IR (reflecting the absorption prop- 
erties of amC dust which strictly decrease with wavelength; cf. 
Fig.D. 

4.3. Photometric amplitudes in different filters 

In Paper II (see Sect. 4.1 there) we outlined that for our mod- 
els dust is the decisive factor for the overall shapes of the SEDs 
and the changing dust absorption also has a strong impact on 
the resulting photometric amplitude during a lightcycle in the 
optical. In contrast, the influence of the changed atmospheric 
structure - being heavily affected by dynamic processes and, 
thus, quite different from the hydrostatic case - on molecular 
features in the resulting spectra is comparably low. For the set of 
models in Table[T]we calculated synthetic lightcurves and com- 
pared the derived amplitudes in Fig.ljQ The dust-free model T 
shows virtually uniform variations (AM ^ 0.6™^^) in all filters. 
All other models exhibit dusty winds and show the characteris- 
tic decrease in amplitudes from the visual towards the IR. This 
gradient steepens for the sequence R — > S — > CI — > C2. This is 
partly a consequence of the increase in Awboi (variations of the 
luminosity input at the inner boundary). However, even more 
important are the corresponding optical depths of the dusty en- 
velopes. The rising amplitudes in Fig.|4] are clearly coupled to 



figure). With the newly emerging dust shell (see also Fig. A.2f 
in Paper II), the optical depths peak at (^boi ~0.6. Since the ab- 
sorption efficiency of amorphous carbon dust has a wavelength 
dependence of approximately oc (see insert of Fig.[TJ, the ef- 
fects are more pronounced in the visual than in the NIR. This 



^ Not consistent with the original dynamical modelling but useful for 
illustration purposes. 

Note that the amplitudes in Fig.|4]are computed from the maximum 
magnitude out of different cycles minus the corresponding minimum 
value and can therefore be larger than amplitudes estimated from the 
two extreme phases of the same cycle shown in Fig.|2] (cycle-to-cycle 
variations). 
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Fig. 4. Amplitudes of the photometric variations for each model 
studied (Table[T]l plotted against the central wavelengths of the 
various filters. 



the more and more pronounced dust formation along the above- 
mentioned model sequence. This leads to a general increase in 
Tdust (Table[TJ, as well as to stronger variations in dust absorption 
during the light cycle (Sect. l4.2T i. 

5. Testing the models against observations of 
C-type l\/liras 

We calculated similar time series of synthetic spectra and pho- 
tometry as in Paper II for the set of models in Table[T] However, 
we restrict the discussion of the modelling results to mean mag- 
nitudes and mean colour indices (i.e., averaged over phase 0boi 
and over different cycles) for the comparison with observational 
data of the corresponding Carbon MirasJ^ Also in the current pa- 
per we used RU Vir as a reference object for model S and over- 
plot the corresponding data in some of the figures below. 

5. 1 . Photometric amplitudes in different filters 

To compare these modelling results with observational findings 
we compiled the data shown in Fig.|5] This plot includes a small 
selection of Carbon Miras from the W06 sample, spanning a rep- 
resentative range in mass-loss rate. While the NIR amplitudes 
were adopted directly from Whitelock et al. (i2006j . we had to es- 
timate the values in the optical range from data in the AAVSC0 
database. Such an estimate can be hampered by the irregular- 
ities in the lightcurves of observed LPVs. For the characteris- 
tic visual amplitudes plotted in Fig.|5] we neglected occurring 
long-term variations (as apparent, e.g. for RUVir in Fig. 2b of 
Mattel 119971 ) that can lead to even higher values (Awvis ~ 8™''^ 
for RUVir). In general, no visual data is available for the ob- 
jects with large M, the amplitude AR of CW Leo was taken from 
the GCVS catalogue (Samus et al. l2012l l. Nevertheless, an over- 
all trend of photometric amplitudes decreasing with wavelength, 
along with a general increase for more evolved stars (more pro- 
nounced pulsation, higher M), can be recognised. 



' Note that a complete set of lightcurves for all standard broad-band 
filters of the Johnson-Cousins-Glass system - as the one for RUVir 
presented in Fig. 7 of Paper II - is not available for other C-type Miras. 

^ phttp : //www ■ aavso . org/data/lcg i 
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Fig. 5. Amplitudes of photometric variations plotted against the 
central filter wavelengths for a selection of targets monitored by 
Whitelock et al. (2006) in the NIR. The estimates for the mass- 
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Wong et al. (l2004l ). and Schoier & Olofsson (|200TT ). 
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dashed line represents the fit given in Eq. (2) of Whitelock et 
al. (2006) for the de-reddened Mira sample. Overplotted are the 
synthetic mean colours of all dynamical models for comparison. 



5.2. Colour indices 

Photometric colours are among the most common properties to 
characterise stellar objects and their overall energy distribution. 
We start the comparison with the colour-colour diagram based 
on NIR filter magnitudes shown in Fig.|6] The observed C-rich 
Miras plotted there cover a large interval in colour, with {J-K)o 
ranging from ^ l.S'^'s to *6.5™s. Whitelock et al. (2006) Hst a 
number of even redder stars with (H-K)o very close to a;4™^S; 
however, such objects lack data in the J-band and cannot be in- 
cluded in Fig.|6] The second reddest object of the well-defined 
sequence in this plot is the well-known dust-enshrouded carbon 
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dynamical models. For the sake of completeness we mark the lo- 
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does not develop a wind (i.e. M = 0). In addition, the empirical 
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Fig. 8. Amplitudes of the /T-band photometric variations of 
the C-type Miras listed in Table 3 of Whitelock et al. (I2006I I 
plotted against the corresponding mean colour indices which 
are corrected for interstellar reddening and given in Table 6 of 
Whitelock et al. (2006). The dotted line was fitted to guide the 
eye. Overplotted are the values resulting from the different atmo- 
spheric models. Note, that model S exactly reproduces the data 
of RU Vir which is, therefore, plotted with a different symbol 
than in other figures. 



star CW Leo (or IRC -1-102 16). To be found at the blue end of the 
sequence (but not shown here; see Fig. 15 in Paper II) is a num- 
ber of C-rich LPVs from the W06 sample with less pronounced 
variability (i.e., most likely SRVs). 

Figure|6] proves that we are able to trace the whole obser- 
vational sequence of Miras by applying the different dynami- 
cal models. Similar to the hydrostatic COMARCS models, all 
the corresponding initial models would be located in a quite 
restricted area in the lower left corner at (H- K)o ~0.5™^^ and 
(7 -//)() Si 1 (not shown here but demonstrated for model S in 
Fig. 1 6 of Paper II), although the stellar parameters of the models 
(Table[T]) differ considerably. Thus, the dusty winds of different 
intensities are the main reason for the C-type Miras stretching 
across the wider colour range in Fig. |6] with redder colour in- 
dices clearly related to higher mass-loss rates. 

5.3. Mass-loss rates 

In general, circumstellar dust grains in the outflows of evolved 
red giants absorb the stellar light in the visual and re-emit it at IR 
wavelengths, as demonstrated e.g. in Fig.[T] On the basis of this, 
(dust-) mass-loss rates of such objects can be derived by fitting 
the redistributed SEDs with the help of dust radiative transfer 
models (e.g. Groenewegen et al. '2009', Riebel et al. 120121) . For 
a review of this topic, a list of modelling approaches, and also 
the related complications, we refer to van Loon (2007 i. For a 
high-quality fit with synthetic SEDs, a sampling of the observed 
SED with a reasonably large number of representative filters is 
desirable. This is often not possible, while it is easier to obtain 
two-colour photometry for the object. Therefore, several stud- 
ies have tried to establish empirical relations between mass-loss 
rates M and various colour indices. A comprehensive list of such 



investigations is given by van Loon (2007 ), also discussing the 
limited accuracy of this method for estimating mass-loss rates. 

Nevertheless, comparing models against observational re- 
sults (e.g. Le Bertre 19971 Le Bertre & Winters 1998 1 in a 
M-colour diagram represents an important test of whether the 
models are adequately describing mass-losing AGB stars. Thus, 
we carried out such a comparison, as was also done for the Berlin 
modelfl simulating dust-driven winds of LPVs with C-rich 
(Le Bertre & Wint ers |1998. Winters et al. I20QO|, as well as 
O-rich (Jeong et al. I2003I I atmospheric chemistry. In Fig.|7]we 
plot the mass-loss rates as derived by Whitelock et al. (2006) for 
their sample of C-type Miras via the mass-loss formula given in 
Jura (119871 ). showing the characteristic rise in NIR colours for 
increasing M. It is interesting that the overplotted empirical re- 
lation found for field carbon stars observed by Le Bertre (1199 7t 
dotted line) lies below the individual data over the whole range. 
On the other hand, the more recent relation for C-rich LMC 
stars adopted from Gullieuszik et al. (120121 dashed line) fits the 
galactic W06 sample better, although these are two independent 
studies. Anyway, the different dynamic model atmospheres, also 
shown in Fig.|7j resemble in general the observational sequence 
of C-type Miras reasonably well. 

5.4. Photometric variations 

The effect of circumstellar reddening due to dusty envelopes de- 
scribed above (Fig.|7]) can in principle also be reproduced quanti- 

' The main difference between their approach and the one used here 
is that we apply non-grey radiative transfer (including molecular opac- 
ities) when computing the dynamic model which has a decisive influ- 
ence on the radial T-p-stmcture. This makes our models suitable to de- 
scribe spectra and photometry of low-mass-loss objects, in contrast to 
the Berlin models. 
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tatively by dust radiative transfer codes, such as DUSTY (Ivezic 
& Elitzur 119971 Ivezic et al. 119991 1. the Groenewegen code 
(Groenewegen [1995 2006 ), 2-Dust (Ueta & Meixner 2003ll, 
or MODUST (Bouwman et al. 2000). Based on simplifying as- 
sumptions for the circumstellar envelopes (concerning the dust 
condensation and composition, the velocity law, the radial den- 
sity structure, the photospheric contribution to the spectra, etc.), 
these codes compute synthetic SEDs for given sets of input pa- 
rameters and prescribed radial structures. Although this may 
provide reasonable relations between colour indices and mass- 
loss rates, another important aspect can by definition not be 
reproduced by such an approach, namely the temporal photo- 
metric variations. In Fig.|8] showing photometric NIR ampli- 
tudes against colour indices, all these models calculated with 
the help of dust radiative transfer codes would line up on the 
abscissa. In contrast, our models have the advantage of account- 
ing for this effect caused by the pulsating stellar interior First, 
they cover the dust envelope and the deep photosphere, pro- 
viding consistent spectra including features due to molecules 
and dust (Fig.[I]l. Second, the pulsation triggered in layers be- 
low the stellar photosphere is simulated by the variable inner 
boundary condition (for details cf. Nowotny et al. 120101 1201 lal l. 
The latter enables us to calculate synthetic lightcurves and over- 
plot our modelling results in Fig.|8] This plot shows the known 
correlation (Olofsson et al. 1 1993b between photometric varia- 
tions and circumstellar reddening (i.e., mass-loss rates). The dy- 
namic model atmospheres are able to reproduce objects across 
the whole range covered by the observational data. 

5.5. Bolometric corrections 

In Fig.|9] we tested the atmospheric models in a different way, 
namely by looking at bolometric corrections (BC) found by ob- 
servational studies. While there is a body of literature available, 
we restricted the comparison in the way that (i) only values 
specifically derived for carbon stars were adopted, and (ii) only 
data given for the standard JK filters were used. 

The first data set plotted in Fig. |9] are the carbon stars from 
the sample of Mendoza & Johnson (,1965; their Tables 1 and 2). 
An empirical fit to these values (not shown in Fig.|9]l was derived 
by Costa & Frogel (fT996] l. given in the CIT/CTIO photomet- 
ric system used by these authors, though. As outlined in Paper I 
(Sect. 4. 1.4), the hydrostatic COMARCS models - a represen- 
tative sample of which is also shown in the figure - are able to 
describe the bulk of the Mendoza objects which are not signif- 
icantly affected by dynamic effects. The mean relation of bolo- 
metric corrections given by Bergeat et al. (2002 , their Table A. 1) 
fits the Mendoza objects rather well, even for colours beyond 
(J-K)q 1.6 where the hydrostatic model approach is not ade- 
quate anymore^ 

The relation provided by Kerschbaum et al. (2010; their 
Eq. (1)), which rests on blackbody-fits to photometric data, ex- 
tends to even larger colour indices. As also noted by Kerschbaum 
et al. (see the discussion in their Sect. 3.3), there seems to be 
a general difference to the bolometric corrections derived by 
Whitelock et al. (2006) for galactic C-type Miras. Overplotted 
also in Fig.|9] the values of the latter authors appear to be most 
suitable for very red C-rich AGB stars which are seriously influ- 
enced by dynamic effects (note that the C-rich LMC objects in 



^ The observed targets of Mendoza & Johnson ( |1965t and Bergeat et 
al. I l2002t which have even bluer {J — K)q colours than the COMARCS 
models of Paper I may be C-rich objects which are not classical AGB 
stars. 
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Fig. 9. Bolometric corrections BCk for carbon stars vs. {J - 
K) colours as compiled from different sources: (i) individual 
measurements for C stars from Mendoza & Johnson (1963] 
squares), (ii) the mean relation for C-type giants of Bergeat et 
al. (2002, full line), (iii) synthetic data based on a representative 
sub-grid of the hydrostatic COMARCS model atmospheres of 
Paper I (crosses), (iv) the fit for nearby field C stars provided by 
Kerschbaum et al. (2010; dash-dotted line), (v) the fit for C star 
candidates in the LMC derived by Riebel et al. (2012 ), (vi) indi- 
vidual measurements for field C-type Miras (grey filled circles, 
blue diamond for RU Vir) taken from Whitelock et al. (120061 
Table 6) as well as the corresponding fit given in their Eq. (10) 
and Table 5 (dashed line). Overplotted are the results for the dy- 
namic model atmospheres of this work. The labels in the upper 
part mark approximate mass-loss rates M [Moyr"'] for the C- 
type Miras of Whitelock et al. (120061 ) as related to a given (7- 
K)q colour in Fig.|7] 



Fig. 12 of Groenewegen et al. 120091 show a similar distribution). 
The models presented in this work are constructed to simulate 
such objects and the corresponding mean values (i.e., averaged 
over phase (/>boi and over different cycles) can be seen in the plot, 
too. 

In observational studies the estimation of bolometric mag- 
nitudes mboi (and, consequently, BC) via integration over pho- 
tometry in various filters is limited by the available data (which 
may be the reason behind the above-mentioned differences be- 
tween different studies, clearly recognisable around (J-K)o ~ 2 
in Fig.|9])- On the theoretical side, we have the advantage that the 
luminosity is an input parameter and bolometric corrections can 
easily be computed as soon as we have the synthetic photom- 
etry. Figure|9] demonstrates that the dynamical models are able 
to trace the sequence of C-rich objects with increasing mass- 
loss rates reasonably well. Some differences at the very red end 
(model C2) may exist; however, the number of observed objects 
around {J-K)() ^6 is also limited. In combination with the hy- 
drostatic models the parabolic trend of BCk vs. (J-K)o for the 
variety of carbon stars can be followed. 
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(m-M)LMC of 18.5"""^. Overp lotted are the C-type Miras adopted 
from Whitelock et al. C2006, cf. their Tables 3 and 6) as well as 
the sequence of dynamic model atmospheres presented in this 
work. 



5.6. Period-luminosity relations 

A different aspect is examined in Fig.[TOl namely the relation be- 
tween the absolute K magnitudes of evolved red giants and the 
periods of their lightcurves. It has been known from observa- 
tional studies of individual objects for quite some time that Mira 
variables follow a period-luminosity relation (PLR; e.g. Feast et 
al. 1989 ). Starting with the seminal paper by Wood (2000 l sev- 
eral authors made use of the data sets provided by extended mi- 
crolensing surveys (e.g. MACHO, OGLE) with the fundamen- 
tal outcome of LPVs constituting parallel sequences in diagrams 
as Fig. [To] These investigations were later extended to different 
stellar systems (e.g. Lorenz et al. 1201 It and other wavelength 
ranges (e.g. Glass et al. |2009l l. 

For illustration we adopted in Fig.[TO] the data for LMC 
variables provided by Ita et al. (2004). The different sequences 
which they could identify empiricially are overplotted and la- 
belled in the plot in the same way as given by these au- 
thors. The prominent Mira variables, pulsating in the fundamen- 
tal mode, are located along sequence C in this plot. Similar 
K-log(P) relations derived by other authors (e.g. Feast et al. 
[19891 Groenewegen & Whitelock [19961 Whitelock et al. [2008] l 
are - within the precision needed for the following comparison 
with our modelling results - quite comparable to the fit plotted 
with a full line. 

It is known that carbon-rich Miras preferentially populate the 
brighter end of sequence C (e.g. Wood ,2000. Ita et al. 2004) . This 
behaviour can be recognised for the C-type Miras of the W06 
sample in Fig.[TOl at least for the more moderate ones (shorter 
periods, bluer colours, brighter in K, lower MLRs). Also the 
models with no or intermediate mass loss (T, R, S) fall right 
onto sequence C constituted by Mira variables with moderate 



outflows. This is not surprising because the stellar parameters 
(Table[T]l where chosen to follow known PLRs, and the mean 
Mfc is related to for these objects. 

However, there is another eye-catching effect in Fig.[TOl 
which becomes apparent for the more evolved objects. As de- 
scribed by Lattanzio & Wood (2004; their Fig. 2.50) and studied 
later in detail by Ita & Matsunaga (201 1 ; their Fig. 2) on the basis 
of observational data, C-rich Miras with pronounced mass-loss 
rates show significant deviations from the practically linear rela- 
tion for fundamental mode pulsators. This effect is wavelength 
dependent, for the K-band such objects appear clearly fainter 
(up to 3™^^) due to the extinction caused by circumstellar dusty 
envelopes (Sect. l4.11 i. The severely reddened targets of the W06 
sample in Fig.[TO|are dimmed by more than 4""^ in K and show a 
distinctive turnoff from sequence C. Figure[TO|illustrates that the 
applied model atmospheres are, in principle, able to reproduce 
this characteristic behaviour found by observations of C-type 
Miras. The models resembling more evolved LPVs with longer 
periods and increasing mass loss (CI, C2) deviate from the lin- 
ear sequence C, although they do not become as faint in Mk as 
the most extreme cases of the observed objects. 

6. Discussion 

6. 1 . The models in the context of a stellar population 

To illustrate how the models generally fit into the global pic- 
ture of a (galactic) stellar population we extracted photomet- 
ric data from the 2MASS survey for Baade's window (Baade 
119631 ). which is known to be one of the low-extinction windows 
towards the centre of the Milky Way (e.g. Dutra et al. ,2002 1) . 
The result of this data mining (cf. App.[Alfor details) is shown 
in the colour-magnitude diagram (CMD) of Fig.[TT] Since the 
Galactic bulge is dominated by an old, red population, a promi- 
nent giant branch can be recognised. It is well populated up to 
the tip of the red giant branch, which is predicted by theory to 
be located at Mboi~-3.9™s or L*29OOL0 (e.g. Salaris et al. 
2002, Cassisi 2012). For orientation purposes the {J-K) colours 
for three representative O-rich objects are indicated at the top 
of Fig.[TT| too, based on the calibrated colour-temperature re- 
lation of solar-metallicity models for late-type giants presented 
by Houdashelt et al. (2000). The sample of Galactic bulge ob- 
ject^ does not contain the kind of stars our models aim to re- 
produce, namely luminous mass-losing C-type Miras. Therefore, 
we overplotted in Fig.[TT]the C-rich field Miras of Whitelock et 
al. (2006) to complete the overall picture. The latter are among 
the most luminous stars in this diagram. The plot illustrates two 
points. First, the luminosities of the models (input parameter, 
cf. Table[r|) were chosen close to reality. And second, the models 
show the same increase in (J-K) compared with the sequence 
of observed C-type Miras. Since this is determined by increased 



' The old but metal-rich stellar population constituting the Galactic 
bulge is dominated by stars fainter than the RGB tip, while only a lim- 
ited number of luminous AGB stars can be found. One would not expect 
to find carbon stars because of (i) the fact of a small intermediate-age 
population, in general, and (ii) the relatively high metallicity. A lack of 
bulge C stars may also be caused by an enhanced oxygen abundance 
(Feast l2007l . Nevertheless, there were studies which claimed to have 
identified C stars in the bulge (e.g. Azzopardi et al. 1991i. Appearing 
quite faint and relatively blue their existence is still puzzling (Whitelock 
et al. 119991 Wahlin et al. 2007 ) and the explanations brought forward 
range from binary evolution (Whitelock ,1999 ) to membership of the 
Sagittarius dwarf galaxy (Ng 11998] 11999) . In any case, this group of 
supposed bulge C stars is consisting of not more than 34 objects. 
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Fig. 11. Colour-magnitude diagram containing objects in 
Baade's window, the bolometric magnitudes were computed 
from 2MASS photometry by using the BCk given in Montegriffo 
et al. (I1998I I. The labels at the top denote effective tempera- 
tures and spectral types of late-type giants for a given cali- 
brated colour as listed by Houdashelt et al. (120001 their Table 4). 
Overplotted are the corresponding data for the galactic C-type 
Miras of Whitelock et al. (2006) adopted from their Table 6, as 
well as the results for all atmospheric models investigated here. 
The vertical dashed line at{J-K)-l.6 marks approximately the 
border between the hydrostatic and the dynamic regime for C 
stars (as discussed in the text). 



circumstellar reddening, this means that the arising mass-loss 
rates are also in a realistic range. 

Frequently the members of a stellar population are divided 
into M- and C-type objects according to a simple colour crite- 
rion. Assuming that all stars redder than a certain NIR colour 
will likely have carbon-rich atmospheric chemistry (Marigo et 
al. |2003j . C stars are singled out photometrically (e.g. Hughes 
& Wo od [19901 Wood|2000 Cioni & Habing 2003, Cioni et al. 
I2OO6I 1. A typical value applied in studies of stellar systems is 
(J-K)- 1.4. However, this has to be regarded as a relatively 
good indication (e.g. Riebel et al. 2012 ) rather than a strict divi- 
sion in chemistry. On the one hand (comparatively warm) C-rich 
objects bluer than this supposed borderline were found by ob- 
servational studies (e.g. the C stars of Mendoza & Johnson[T965] 
in Fig.|9]l as well as in our hydrostatic model grid presented in 
Paper I (all models there have (7-7*^) < 1.6'™°). On the other 
hand there are evolved O-rich objects known which are redder, 
such as the dust-enshrouded M-type Mir a IRC- 20197 (IW Hya), 
which exhibits according to Le Bertre (1993, see their Fig. 3) 
a (J-K) of 3™° and more. Other examples can be recognised 
among the Galactic bulge Miras of Groenewegen & Blommaert 
(12005) plotted in Fig. IB. II or in studies of very evolved OH/IR 
stars (e.g. Marshall et al. 12004) . What can be concluded from 
our modelling, though, is that (J-K) ^1.6 seems to be the start- 
ing point from where the dynamic aspects - pulsation and es- 
pecially the dust effects - become relevant for C stars (marked 
with a dashed line in Fig.fTTTi. All of the hydrostatic dust-free 
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Fig. 12. Colour-magnitude diagram of the Galactic bulge ob- 
jects together with the reddening corrected mean magnitudes 
and colours of field Carbon Miras from Whitelock et al. (2006] 
Table 6) shifted to an absolute scale Mfc according to the dis- 
tances given by the same authors. Overplotted are the mean 
colours and magnitudes of the series of dynamical models. 
Marked are, in addition, the magnitude of the RGB tip esti- 
mated by using the empirical relations of Valenti et al. (120041 
Eq. (9)) under the assumption of [Fe/H] ^ for the Galactic 
bulge (Zoccali et al. 120031 ). as well as the magnitude of the be- 
ginning of the TP- AGB for a solar-like star as given by Marigo et 
al. (,2008) . The labels in the upper part mark approximate mass- 
loss rates M [Mq yr"'] for the C-type Miras of Whitelock et al. 
(120061 ) as related to a given {J-K)o colour in Fig.|7] 



COMARCS models of Paper I have bluer colours. One of the 
brightest and nearest optical carbon stars is the extensively stud- 
ied TX Psc which only shows small photometric variability and 
a weak mass loss (e.g. Gautschy-Loidl et al. 120041 ). Its 2MASS 
photometry (Cutri et al. 120031 ) transformed using the relations 
of Carpenter (2001) results in (J-K)Kl.n, placing TXPsc 
only slightly beyond the assumed border. Clearly redder than 
{J-K) = l.6 we find in Fig. [TT] the Mira variables with pro- 
nounced variability and significant stellar winds, both from the 
observational sample of Whitelock et al. (120061 ) as well as the 
corresponding models of our study. Thus, a selection of red 
stars based on this colour criterion will essentially capture non- 
hydrostatic, reddened C stars. 

Figure[T2] shows a similar colour-magnitude diagram for 
our combined galactic population - including the bulge objects 
in Baade's window and the C-rich field Miras adopted from 
Whitelock et al. ( 120061) - with the dynamic model atmospheres 
overplotted. In contrast to the previous Fig.[TT] where the car- 
bon Miras show luminosities steadily increasing with (J-K)o, 
the same objects show here a characteristic dimming in Mk due 
to the shift of the SED maximum (cf. Sect. 14. 11 ). reflecting the 
growing importance of circumstellar dust. Similar CMDs were 
documented by Zijlstra et al. (120061 ) and Lagadec et al. (120071 
their Fig. 2) for the Magellanic Clouds and the C stars observed 
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in these systems with the Spitzer Space Telescope. It can be seen 
that the models are able to trace the sequence of Carbon Miras 
from slightly above the RGB tip along the curved branch towards 
higher {J-K)o colours. This sequence is related to increasing 
mass-loss rates as denoted in the figure{3 There may be some 
deviations in the amount of dimming in Mk when model C2 is 
compared to the very reddened observed objects. On the other 
hand, the observational sequence is also only sparsely populated 
at {J-K)o x!6'™^. Further work concerning these discrepancies 
is required to investigate reasons and remedies. 

6.2. (C-0) versus [Fe/H] 

For the comparison of our modelling results with observational 
findings we also adopted empirical relations that were origi- 
nally derived for C-rich objects in the Large Magellanic Clouds 
(LMC). The first case can be found in Fig.|7] where the mod- 
els were contrasted with the objects of the W06 sample in the 
M-{J-K)-plane. Also shown in the figure is the LMC relation 
derived by Gullieuszik et al. (120121) . which fits both kinds of 
data quite well. The second case concerns the bolometric correc- 
tions compiled in Fig.|9] It appears that the fit derived recently by 
Riebel et al. (2012) for candidate C-rich AGB stars in the LMC, 
which is overplotted in Fig.|9] provides a relatively good over- 
all representation of the bolometric corrections BCk combined 
from various data sets (with differences recognisable for the very 
red objects). 

The agreement between, on the one hand, our models with 
solar-like composition and the corresponding galactic observa- 
tional data (esp. the field C-type Miras of Whitelock et al. 120061 ) 
and, on the other hand, the LMC relations is noteworthy as the 
metallicities of the Milky Way and the Magellanic Clouds differ 
(general values of [Fe/H] a; -0.5 and lower were estimated for 
the LMC). This suggests that a difference in metallicity is not 
the decisive factor for mass loss and, thus, the reddening of the 
objects. In fact, the amount of excess carbon (C-O) - meaning 
available C atoms that are not locked in CO molecules - appears 
to be more relevant to the dust formation, the properties of the 
developing wind and the resulting circumstellar reddening. This 
supports the theoretical results of Mattsson et al. (I2008t . 

6.3. Extremely reddened C stars 

Frequently, studies dealing with photometric data obtained for 
the Magellanic Clouds by the Spitzer SAGE survey (e.g. Vijh 
et al. 120091 Boyer et al. 12011) follow the scheme introduced 
by Blum et al. (120061) to divide the found AGB objects into 
three distinct classes tagged "O-rich", "C-rich", and "Extreme". 
Recently, efforts were made to interpret the SEDs of these ob- 
jects (defined by photometry compiled from various catalogues; 
cf. Srinivasan et al. 120091 ) with the help of a model grid (Sargent 
et al. l20m Srinivasan et aL lSOTTI Riebel et al. 2012^. Thi s SEP 
fitting led to the conclusion (cf. Fig. 1 of Riebel et al. 120121 ) 
that the vast majority of the identified "extreme AGB stars" are 
obscured C-rich objects with high mass-loss rates, while only 
very few of them might be bright O-rich objects undergoing 
heavy mass loss (OH/IR stars). In this light, the border between 
"C-rich" an d "ext reme" AGB stars (set to (7-^:) ss2.1™''s in 
Blum et al. 120061) appears to be somewhat arbitrary. The lat- 



For orientation purposes, a selection of observed targets is also 
marked in Fig.[T2] with mass-loss rates as derived by Whitelock et al. 
([2006}: V Oph (M« 1.1 x 10"^ Meyr"'), RU Vir (2x10-"), FX Ser 
(7.6 X 10^"), RW LMi (1.1 x 10^'), and CW Leo (2.2 x 10^^). 



ter objects should not be regarded as a distinct class of targets 
(as suggested by the classification scheme) but rather as C-stars 
with severe reddening due to circumstellar dust. From the re- 
sults presented in this work we would assume a natural transition 
from "C-rich" to "extreme" simply caused by increasing mass- 
loss rates. Such IR-C-stars with colours as red as (J-K)^6 are 
also known from previous photometric investigations as, for ex- 
ample, the 2MASS survey of the LMC (Nikolaev & Weinberg 
|20U0] ) or the Galactic bulge (Cole & Weinberg [2002] ). 



7. Summary and conclusions 

The aim of this work was to reproduce the observable variety 
of mass-losing galactic Miras of carbon-rich atmospheric chem- 
istry with the help of sophisticated model atmospheres. In partic- 
ular, we wanted to see if our modelling approach is able to sim- 
ulate the photometric properties of the comprehensive sample 
of well-characterised C-type field Miras observed by Whitelock 
et al. (120061 ). These objects form - more or less pronounced - 
sequences in various kinds of diagrams. Several properties char- 
acterising the stars (e.g. luminosities, colours, pulsation periods, 
photometric amplitudes) show on average an increasing trend 
along the sequences. This could suggest an evolution of C-Miras 
along this series of observed targets. However, one should keep 
in mind that the W06 sample of field Miras might contain a mix- 
ture of different masses and ages, making the interpretation as 
an direct evolutionary sequence arguable. Nevertheless, a gen- 
eral trend from less to more evolved objects is recognisable, with 
the increasing mass loss being one of the most significant effects 
causing circumstellar reddening. 

For our photometric modelling we used time-dependent dy- 
namic model atmospheres that simulate the dynamic photo- 
spheres, as well as the pulsation-enhanced dust-driven winds 
of C-rich LPVs (Hofner et al. 2003 Mattsson et al. 120101 ). It 
should be noted that this self-consistent approach differs from 
frequently used dust radiative transfer codes to model dust en- 
velopes of mass-losing AGB stars. One major difference is the 
fact that the development of a stellar wind is the result of the 
chosen set of fundamental parameters, i.e., the mass-loss rate M 
is not an input parameter but actually an outcome of the mod- 
elling. Another advantage is that we are also able to compute 
variations in spectroscopic and photometric properties, which is 
crucial for intrinsically variable objects such as Miras. 

A set of five models with quite different parameters was used 
to represent the diversity of Carbon Miras found in observational 
studies. Based on the radial atmospheric and wind structures we 
computed (consistent with the preceding modelling) synthetic 
spectra and photometry for standard broad-band filters. In a first 
step we presented low-resolution spectra for different models, 
illustrating that our modelling method is able to describe a di- 
versity of objects ranging from only moderately pulsating LPVs 
without any dusty outflow to evolved Miras with pronounced 
variability and optically thick circumstellar dust envelopes. In 
this way, our modelling approach bridges the gap between (i) the 
Australia-Heidelberg-type models for dust-free pulsating atmo- 
spheres without wind, and (ii) the Berlin-type models for dust- 
driven winds of highly dust-enshrouded objects (see Paper II and 
Sect. 3 of Nowotny et al. 12005a] for a more detailed description 
and references). Our models can reproduce spectra of C-rich 
AGB stars, for which the respective contributions by the pho- 
tosphere and the dusty envelope are of varying relevance. This 
also includes objects with significant mass loss but no optically 
thick winds. 
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The applied models compare well to the observational data 
and can (quantitatively and simultaneously) reproduce the ob- 
served values and ranges for absolute magnitudes, colour in- 
dices, photometric variations, mass-loss rates, and bolometric 
corrections. We were able to retrace the sequence of galactic 
C-type Miras from the "blue end" (stars slightly beyond the tip 
of the RGB, modest variations, no outflow) to the "red end" 
(evolved red giants, large photometric amplitudes, pronounced 
circumstellar reddening due to heavy mass loss) in diverse di- 
agrams: CMDs, two-colour diagrams, M vs. (J-K), AK vs. 
(J-K), or K-log{P) diagrams. From this we can conclude that 
(i) the input parameters of the models were chosen realistically 
(e.g. Li, + Aup + /l — » AK reproduced), and (ii) the stellar winds 
emerging in the modelling scenario have properties resembling 
weU the corresponding observed objects (e.g. M + u^xp+fc — > 
red (J-K) colours reproduced), suggesting an adequate mod- 
elling procedure. The combination of both enables us to repro- 
duce the diversity of observed populations of C-rich Miras. 
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Appendix A: Various comparative data 

In Sect.|6?T](and App.|B]l we use phot ometry of the 2MASS All- 
Sky Point Source Catalog (Cutri et al. 2003) for a section of the 
Galactic bulge, which is available from the NAS A/IPAC Infrared 
Science Archive (IRSA)El Photometric data were downloaded 
for objects of a field located in Baade's window, we extracted 
photometry for an area of Ao- y. AS ^ 0.5°x 1°. Only stars ful- 
filling cTjHK < 0.05™^^ were chosen, the correction for interstel- 
lar reddening was carried out according to the extinction values 
given by Dutra et al. (|^002 ). The 2MASS photometry had to be 
transformed to our Bessell system using the equations given in 
Carpenter (120011 see their App. A). To convert the /^-magnitudes 
of 2MASS to bolometric magnitudes (Fig.fTTTl we made use of 
the bolometric corrections BCk as a function of (J-K) derived 
by MontegrifFo et al. ( I1998I I. To this end we fitted the values tab- 
ulated in Montegriffo's Table 3, which are labelled 'metal-rich 
stars' . The resulting regressiorQ was then applied to the 2MASS 
data. 

In addition, we used a sample of O-rich Miras in the Galactic 
bulge identified by Groenewegen & Blommaert (2005) while 
investigating the lightcurves produced by the OGLE-II survey. 
The authors cross-correlated their objects with the 2MASS cata- 
logue leading to a subsample of 1619 LPVs which are plotted in 
Fig. IB. II Dereddening was done following the values for inter- 
stellar extinction for each OGLE-II field as listed in Matsunaga 
et al. (I2005I I. For the conversion between 2MASS and Bessell 
filters we applied again the relations of Carpenter (1 2001 ). 

The reddening-corrected photometric data (apparent magni- 
tudes) for all bulge objects were shifted to absolute magnitudes 
by assuming a distance modulus of (m-M) - 14.7™"° for the 
centre of our galaxy (e.g. VanhoUebeke et al. 120091 1. 

For the comparison in Fig. [10] we made use of the data for 
variable stars in the Large Magellanic Cloud (LMC) presented 
in Ita et al. (2004). The authors analysed OGLE lightcurves and 
provide estimates for the periods for a large number of objects, 
as well as the corresponding JHK photometry. We dereddened 
the data with the value for Ak given by Feast et al. ( 11989b and 
transformed the near-IR photometry from the LCO to the Bessell 
system with the help of Carpenter (2001 1. We also adopted the 
P-K relations fitted by Ita et al. ( 2004 ; cf. their Table 3) and over- 
plotted them in Fig. [10] As a distance modulus for the LMC we 
assum ed (m -M) - 18. S™"*" (e .g. Gib son 12000. Freedman et al. 
I200n Alves 120041 Laney et al. l20T2l ). 

The relation between mass-loss rates M and (J-K) colours 
empirically found by Le Bertre (1997) and plotted in Fig.]7]is 
given in the ESO photometric system and was transformed to 
the Johnson-Glass system according to the equations given in 
Bessell &Brett(1988). The other relation given in Gullieuszik et 
al. (2012 i their Eq.(2)) was transformed from the 2MASS system 
by applying the conversion of Carpenter (2001). 

The magnitudes of the solar-like hydrostatic COMARCS 
models - a sub-grid is used in Fig. [9]- are of the same kind as 
the corresponding ones of the dynamic models (Paper I). 

Finally, we compiled various data sets from different sources 
for the comparison of bolometric corrections (BCs) shown in 
Fig. [9] The observational results of Mendoza & Johnson (1965 ) 
are akeady given in the Johnson system, no conversion is needed 
therefore. The fit to the values of an extended sample of C-type 
giants provided by Bergeat et al. (120021 their Table A. 1) is based 
on observed data collected by the authors from the literature. 



" http : //irsa . ipac . caltech . edu/ 

'2 BCk = O.m + 4.545 ■ (J - K)-2JS6- (J -K)- + 0JS6- (J -Kf 



According to Knapik et al. (119971) these measurements were ob- 
tained in the "Arizona system or close to it" (i.e., Johnson) and 
can thus be compared directly (Bessell & Brett 1988, Table I) 
to the models. The fit to the data of a sample of nearby field 
C stars provided by Kerschbaum et al. (120101 ) in their Eq. (1) 
is based on photometry in the new ESO system, which is very 
close to the Johnson system according to the authors (see also 
Bouchet et al. 1991). The fit for C-rich AGB candidates in the 
LMC given by Riebel et al. ( 120121) in their Table 7 was trans- 
formed from the 2MASS system to our standard system with the 
help of Carpenter ( 2001 ). 



Appendix B: Loops in colour-magnitude diagrams 

In the previous Sects. [5] and [6] we analysed the applied models 
on the basis of mean photometry. Here we want to outline how 
typical C-type Miras vary within CMDs constructed from NIR 
photometry. 

In Fig. lB.ll we show the same combined population of galac- 
tic objects as in Fig.[TT] but with absolute K magnitudes on the 
ordinate instead of bolometric ones. This enables us to overplot 
also the Miras (most likely 0-riclJ2I) found by Groenewegen & 
Blommaert (120051 ) in the 49 OGLE fields across the Galactic 
bulge. For a description of the necessary post-processing of 
the 2MASS photometry we refer to App.|A] The fact that the 
2MASS survey provides 'only' single-epoch data (instead of a 
time series that could be averaged) introduces some scatter in 
for these variable stars. Still it is clear that the latter objects are 
all located on the extension of the pronounced red giant branch 
towards even brighter Mk in Fig. IB. II The C-type Miras of the 
W06 sample form again the noticeable sequence from the bright 
end of the red giant branch towards redder colours with (J-K)o 
up to ~6™"°. Also plotted in Fig. IB. H is the synthetic photome- 
try for model S, illustrating the difference between the dust-free 
hydrostatic initial model and the developed dynamic model in a 
CMD. 

Figure lB72] (showing a close-up of Fig. IB. lb illustrates the 
temporal variations, with model S and the C-type Mira RU Vir 
as examples. Both occupy the same range of the CMD and show 
variations on the order of 1""° in K as well as (J-K\). Apart 
from the recognisable cycle-to-cycle variations the model shows 
for every pulsation period similar broad loops throughout the 
light cycle. It appears bright and blue at maximum light (phase 
0boi ~ 0.0), while it is faint and red around minimum light 
(0boi ~ 0.5). As described in Paper II, we cannot plot such loops 
for RU Vir directly as the observed data points were obtained 
with an insufficient sampling rate over several periods. Instead, 
we merged the measurements from different periods into one 
combined light cycle. Based on the sine fits (one Fourier 
component) to the resulting lightcurves (see Fig. 8 in Paper II) 
we can derive simulated variations of RU Vir, which are shown 
with a dashed line in Fig. IB. 21 The range of this simulated loop 
is narrower than the range covered by the individual data points. 
As noted in Paper II (cf. the comparison of original (J-K) data 
vs. the corresponding simulated colour curve in Fig. 12 there), 
this can be explained by the occuring cycle-to-cycle variations 
in the lightcurves. Nevertheless, the general behaviour of 
looping through the CMD, with extreme phases in the same 
order (0boi ~ 0.0 at the upper left end, 0boi ~ 0.5 at the lower 
right end), is similar to the model. Again we find that the 
model loops counter-clockwise in the CMD of Fig. IB. 21 while 
RUVir follows a clockwise movement. We demonstrated in 
the appendix of Paper II how a (relatively small) phase shift 
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Fig.B.l. Colour-magnitude diagram containing observational 
data compiled from different sources: (i) 2MASS photometry of 
stars in Baade's window (black dots), (ii) 2MASS photometry 
of Mira variables in the Galactic bulge (turquoise dots) iden- 
tified by Groenewegen & Blommaert (I2005I I. and (iii) redden- 
ing corrected mean magnitudes and colours of field C-Miras 
(grey filled circles, blue diamond for the individual RU Vir data) 
from Whitelock et al. (120061 Table 6) shifted to an absolute 
scale Mk according to the distances given by the same authors. 
Overplotted are the modelling results for model S, i.e. the loca- 
tion of the hydrostatic initial model as well as the average of sev- 
eral phases of the different cycles in Fig. IB. 21 The colour-code is 
the same as in Fig. 16 of Paper II. The box (dotted lines) marks 
the range covered in Fig. IB. 21 



between the lightcurves obtained in the different filters can lead 
to a change in the sense of rotation. Interestingly, this behaviour 
has only scarcely been investigated directly by observational 
studies. Payne-Gaposchkin & Whitney ( 119761 their Figs. 5abc) 
studied a large sample of LPVs and found them - averaged 
over certain groups - to vary in CMDs very similarly to our 
results. Unfortunately, they show only loops for stars of spectral 
types M and S (clockwise), while the corresponding data for 
C-rich objects was according to the authors afflicted with large 
scatter. Eggen (!1975', see his Figs. 18, 24, and 30) monitored 
large-amplitude variables (mostly M-type) and plotted their 
variations in the Mboi -(/?-/) -plane. Although the sampling of 
these time series is limited, the stars seem to perform loops 
with the sense of rotation being mostly counter-clockwise. A 
slightly different approach was pursued by Wing (1967), who 
obtained a spectro-photometric time series over » 2.5 years for 
about 25 Miras from which he could derive oxide band strengths 
and temperature estimates. According to Wing (priv. comm.) 
the objects showed quite a variety of behaviours. The results 
for only five of these objects were then published in Spinrad 
& Wing (119691 their Fig. 9), showing clockwise as well as 
counter-clockwise variations. 
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Fig.B.2. Colour-magnitude diagram containing individual ob- 
servations of RU Vir adopted from Table 2 of Whitelock et al. 
(120061) and shifted to an absolute scale Mk by applying the dis- 
tance modulus derived by the same authors (cf. their Table 6). 
The dashed line marks the average variation of RU Vir during the 
light cycle as derived from the Fourier fits in Fig. 8 of Paper II. 
Overplotted are the corresponding synthetic photometric data of 
model S. Three different pulsation cycles are plotted with the 
same symbols, colour-coding, and labels as in Fig. 13 of Paper II. 
The arrows mark the directions of how the objects pass through 
the loops during a light cycle, selected phases are labelled (con- 
vention for 0boi / 0v as in Paper II). 



Note that in Fig. 3 of Nowotny et al. (1201 Ibl) we plotted a 
colour-magnitude diagram Mk vs. (V-K)o, showing the com- 
bined contents of Figs. IB. ll and IB. 21 from here. Although this 
differing filter combination provides a valuable colour index be- 
cause of the broad basis in wavelength, it is generally harder 
to obtain simultaneous measurements in these two filters. For 
most of the objects in the W06 sample there is no informa- 
tion concerning their brightnesses or even light variations in the 
V-band and they are not included in this plot. For a quite lim- 
ited sub-sample of the C-type Miras we could estimate a mean 
visual magnitude from the photometric time series available in 
the AAVSO database, providing a hint on a mean (V-K)(). This 
is only possible for stars that are bright in the visual and not 
strongly influenced by circumstellar reddening. Thus, only ob- 
jects with M<2x lO^^Moyr"' could be studied. The C-rich 
Mira RU Vir is the only target for which simulated variations in 
the MA:-(y -A')o-plane could be derived from the Fourier fits of 
the lightcurves (see Fig. 8 in Paper II). Both RU Vir and model S 
show a relatively similar behaviour as found here (loops), with 
the averaged colour of the latter being significantly redder, 
though. 

Appendix C: Extreme dust formation 



In Sect. 14. 21 we illustrated how the time-dependent dust forma- 
tion in the outer layers of a C-type Mira effects the resulting 
lightcurves of such an object. This was based on a dynamic 
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Fig. C.l. Upper panel: Movement of the outer layers of a spe- 
cific dynamic model atmosphere with dust formation occuring 
every second pulsation period (App.|C)- Colour-coded is the 
degree of condensation of carbon atoms into amorphous car- 
bon dust particles. Lower panel: Synthetic lightcurves for the 
broad-band filters V and K (filled circles), together with the cor- 
responding photometry for which the dust opacities were ne- 
glected (open circles). Note that the K-lightcurves were shifted 
by 3™^ to better fit into one panel. 



model atmosphere with relatively moderate dust formation pro- 
cesses. Here we show a different modeQ adopted from the re- 
cent grid of Mattsson et al. ( 1201 Oi l which represents the quite 
rare group of strictly 'double-periodic' models (cf. the discus- 
sion in Hofner et al. l2003l l. As it can be seen in the upper panel 
of Fig. lC.ll only every second pulsation period a new dust shell 
emerges and propagates outwards due to the radiation pressure 
acting upon the dust grains. Figure lC2l shows that the radial 
structure of the model replicates in the inner, dust-free layers 
from cycle to cycle for a given phase 0boi- However, the layers 
outwards from the dust-forming region at a; 2-3/?* can differ 
significantly because of the extreme dust formation. The emer- 
gence of the quite narrow dust shell is also reflected in the pho- 
tometric variations plotted in the lower panel of Fig. lC.ll (see 
also the discussion in Sect. l4.2b . While the K-lightcurve follows 
approximately the sinusoidally changing luminosity input, the 
model exhibits especially in the visual a characteristic lightcurve 
which is reminiscent of 'Miras with secondary maxima' as dis- 
cussed e.g. by Lebzelter et al. (I2005al l. It would be interesting to 
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Fig. C.2. Radial atmospheric structures for the dynamic model 
used for Fig. lC.ll and discussed in App.lCl Shown are two 
characeristic phases (marked with dotted lines in Fig. IC.lt where 
dust formation takes place (0boi = 2.573) or not (0boi - 1.585). 



make a more detailed comparison with observed targets show- 
ing a similar behaviour, as for example the cases presented in 
Wood (IT9991 object 7308. 11 3 in their Fig. 2), Groenewegen et al. 
(20 041 obj ect 050709.50-685849.4 in their Fig. 9), Lebzelter et 
al. (2005a; R Nor in their Fig. 8), Lebzelter et al. (2005b, LWIO 
in their Fig. 2) or Lebzelter (i2011: , ASAS 201445-4659.0 in their 
Fig. 1), in the future. 
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